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Burgers turbulence supported by white-in-time random lorcing at low wavenumbers is studied 
analytically and by computer simulation. The peak of the probability distribution function (pdf) 
Q{^) of velocity gradient ^ is at ^ = where 5/ is a forcing parameter. It is concluded that 

Q(^) displays four asymptotic regimes at Reynolds number > 1: (A) Q{^) ~ C/^C exp(-^^/3C/) 
for C > C/ (reduction of large positive ^ by stretching); (B) Q{^) ~ C/ICT^ for R^^^O 
(transient inviscid steepening of negative ^); (C) Q(C) ^ for R^^^^f ^ ^ R^f (shoulders 

of mature shocks); (D) very rapid decay of Q for — ^ > 0{R£,f) (interior of mature shocks). The 
typical shock width is 0{l/Rkf). If R~^^^ 3> rkf 2> , the pdf of velocity difference across an 
interval r is found be P{Au,r) oc r~^Q{Au/r) throughout regimes A and B and into the middle of 
C. 



I. INTRODUCTION 

A number of papers in the past several years have 
treated the physics of turbulent solutions of Burgers 
equation [p]-p^. A variety of results on statistics of 
solutions, sometimes contradictory, have been reported. 
Burgers turbulence driven by forcing that varies infinitely 
rapidly in time (white-in-time) was first systematically 
studied by Chekhlov and Yakhot ^j^- The present pa- 
per is concerned with Burgers turbulence supported by 
white-in-time forcing that has a compact wavenumber 
spectrum. 

Two processes act upon the velocity field injected by 
a forcing term: the self-advection of the velocity field 
steepens negative velocity gradients and reduces positive 
gradients. The viscous term relaxes the curvature of the 
velocity field. The effects of each term are easy to treat 
in isolation, but the combination of the two terms poses 
nontrivial difficulties. 

At large Reynolds numbers, the self advection tends to 
produce sawtooth structures with smooth ramps of gen- 
tle positive velocity gradient and narrow shocks of strong 
negative gradient. The shock width is 0{iy/Au), where 
Au is the jump in velocity across the shock and v is vis- 
cosity. Shocks can interact. A strong shock can move 
across the domain and swallow weaker structures in its 
path. 

The present paper offers physically motivated approx- 
imations on terms in the exact equations of motion for 
the probability distributions of velocity gradient and ve- 
locity difference and examines their consistency. We as- 
sume that the forcing has compact spectral support at 
low wavenumbers and is white in time. The results are 
tested against computer simulations and the relation to 
other theoretical approaches is discussed. Four asymp- 
totic regions are predicted at large Reynolds number: a 



region of large positive gradient where the gradient prob- 
ability density decreases very rapidly; a region of inter- 
mediate negative gradient where the density follows a —3 
law; a region of larger intermediate gradient where the 
density follows a —1 law; and finally an outer region of 
negative gradient where decay is very rapid. The power- 
law regions are mediated by transient advective steep- 
ening of gradients and the shoulders of mature shocks, 
respectively, while the outer region represents gradients 
within strong shocks. 



II. STATISTICAL EQUATIONS FOR THE 
VELOCITY GRADIENT 



Burgers equation with forcing is 
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where u{x,t) is a one-dimensional velocity field and 
f{x,t) is a forcing term. The left side of ([^) is the La- 
grangian time derivative of u, measured along a fluid- 
element trajectory. We shall assume that the forcing 
is white in time, statistically homogeneous and station- 
ary, with compact spectral support concentrated about 
a wavenumber kf. Let 
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{fx{x,t)fxix, s)) ds, 
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where u{x,t = 0) = and ( ) denotes ensemble average. 
A characteristic forcing strain rate and Reynolds number 
induced by the forcing may be defined by ^/ = B^^^ and 
R = ^//(i/fcp. The steady-state values of rms velocity 

1 /2 

Urms — i}^'^) ^-iid typical shock jump induced by the 
forcing are both 0{£,f/kf). The typical shock width is 
V /urms and the typical velocity gradient S, = Ux within a 
shock is i?^/. 
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Differentiation of (||) yields 

£,t + U^x = + l^£,xx + fx 



(3) 



Tlie term in represents advective gradient intensi- 
fication or diminution along Lagrangian trajectories. An 
exact equation of motion for the probability distribution 
function (pdf) Q{^) of ^ follows from (||): 



de ■ 



(4) 



Here H{£^) = {Cxx\0 denotes the ensemble mean of £,xx 
conditional on fixed ^. This relation was derived in Q for 
zero forcing by following probabilities along Lagrangian 
trajectories. The B term in (Q) expresses in standard 
fashion the outward diffusion of probability due to white- 
in-timc forcing. 

The ter m in (||) plays two opposed roles. If ^ < 0, 
it intensifies the gradient but, at the same time, squeezes 
the fluid and thereby decreases the measure along x as- 
sociated with an interval d^. If ^ > 0, the gradient is 
decreased but measure is increased by stretching of the 
fluid. The intensiflcation or diminution of gradient is 
expressed in (||) by the d{£,^Q)/d£_ term and the rate of 
change of measure is expressed by the term. An iden- 
tity for homogeneous flelds. 
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{Hx\OQ) = CQ, 



(5) 



gives an alternative expression for the rate of change of 
measure. 



III. APPROXIMATIONS FOR ASYMPTOTIC 
RANGES 

The limit R ^ oo suggests the existence of several 
asymptotic ranges of Q(^), whose form can be found if 
physically-based approximations to the terms of (^) are 
valid. We shall outline these approximations compactly 
in this Section, and later address some questions and 
paradoxes that arise. 

First, consider places in the flow, at a given time, where 
forcing has produced an extraordinarily large positive 
value ^ 3> Stretching should very quickly flatten 
such regions, suggesting that the v terms in (|^) and (Q) 
can be neglected in comparison to the other terms in a 
statistically steady state. The Q equation then reduces 
to 



B 



(6) 



The general solution involves exp(— ^■^/3i3). The solution 
that vanishes at ^ = -|-c» is 



(7) 



where C+ is a dimensionless constant and the inequality 
defines the range in which we hope the neglect of viscous 
effects is justified. The exponential factor in (|^) was first 
found by Polyakov but with a different prefactor. 

Next, consider places where ^ < and ^/ <C — ^ <C 
R^f. The first inequality implies that the gradient is 
steepening with a time constant <C 1/^/ and the second 
inequality makes |^| small compared to typical gradients 
found within a shock. We therefore neglect the and B 
terms in (0), yielding in steady-state 
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(8) 



The solution of physical interest is 

Qio^c^seficr' i^f«~^«R'/%), (9) 

where C_3 is another dimensionless constant. Again, 
the inequalities express the range in which we hope 
our approximations are valid. The second inequality 
stated in (^) is stronger than called for by the argu- 
ment just given. This is because the transient steep- 
ening contribution (|^) to Q is overshadowed in the range 
R^^^S,f <^ ^ R£,f by contributions from the shoul- 
ders of the strong, quasi-equilibrium shocks that decay 
with time constants 0(1/^/). 

There is a simple physical explanation for (|9|). Con- 
sider a place where forcing has produced negative ^ with 
—1^ = 0{^f ) over a region of size 0{l/kf ) where £,xx is 
small enough that viscous effects are negligible. Steep- 
ening of the gradient will carry the gradient to larger 
negative values and the time in which it is 0(|^|) will 
be 0(1/|^|). Simultaneously, the measure will shrink to 
Thus this process will contribute 0(1/|C|2) to 
the probability of finding the gradient at 0(|^|). 
This implies Q(0 oc 

An equilibrium single-shock solution of (|l|), (|^) is 
u{x) = — Us tanh(usa:;/2;/). 
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The jump across the shock is 2us. The measure on x of 
places with gradient between ^ and is df/ICa:!. By 

differentiation of (^0|), the contribution of this shock to 
Q{£,) is then 



(11) 



If Us = 0(urras) and |^| <C ^0', ^^^^^ simply Q(^) oc 

Now consider a train of such shocks, spaced 1/fc/ apart, 
with a uniform gradient 2uskf added so that there is no 
secular change of u. It follows that the normalized total 
contribution to Q{^) for ^/ < |^| < R^j is 0(l/i?|C|). 
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The significance of R here is that it measures the ra- 
tio of shock spacing to shock width. For |^| ^ i?^/^^/, 
this contribution to Q is overpowered by (^), but for 
1^1 ^ R^^'^^f, it overpowers (||). Thus the third asymp- 
totic range is 



iR'/% « « R^f) 



(12) 



where C_i is yet another dimensionless constant. This 
result is unchanged in form if the shocks of fixed jump 
2us are replaced by shocks with a distribution of jump 

values Us = OJUrms)- 

The form dl3) may also be obtained by approximate 
dynamic analysis that follows the evolution of those re- 
gions of negative gradient that evolve into the shoulders 
of an equilibrium shock. An example is the mapping clo- 
sure for Burgers shocks [||. 

Equation ( pj] ) may be integrated over a distribution 
Ps{us) of shock jumps to get the form of Q{(,) for 
— ^ > 0{R^f), where contributions from within the 
shocks must dominate Q. Thus, 



Ps{Us)dUs 




(-e > 0{R^f)), 
(13) 



where Cs is a dimensionless constant. Ps(us) oc 
Ms exp(— iaMg/w^^s), where a is a dimensionless con- 
stant, is one of the forms for which the integration in 
( p^ ) is analytic. It yields 



exp(-a|C|i?/0) (-e > 0(i?e/)). (14) 



It happens that the form ( p^ is also obtained by apply- 
ing mapping closure to (jl|) and (||), both in the case of 
free decay and in the present forced case. 

Note that ( p^ ) gives an exponential fall-off while the 
assumed Ps{us) falls off as a Gaussian. This is because 
the gradient in the center of a shock is —u1/2v. The 
shock width v/ug decreases as Ug increases. 

Several authors have concluded that the tail of the pdf 
of shock jump falls off more rapidly than that of the 
Gaussian forcing; in particular, \uPs(us) oc —{ug/urms)'^ 
at large Us [§|,||J|]- The reason is that the decay rate 
of a shock is 0{uskf), which increases with increase of 
Ms. The corresponding prediction for the gradient pdf is 
\nQ{^) cx — (ICI^/C/)'^^^ at large negative 5, instead of 
(0)- 



IV. THE CONDITIONAL MEAN OF 
DISSIPATION 

The conditional mean H{^) embodies all the difficulty 
in solving (El). We wish now to examine whether this 



mean behaves in a way that realizes the asymptotic 
R^ oo ranges of Q{^) proposed in Sec. III. 

The steady-state asymptotic behavior of H{£^) at — ^ » 
R^f can be inferred immediately if the B term can be ne- 
glected in that range. H must balance the remaining two 
terms, which represent steepening and loss of measure: 
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(15) 



Note that statistical homogeneity requires 
/_!lQ(OC(^C = 0. If g(0 falls off faster than alge- 
braically as f — oo, then the term in ( [T^ ) is domi- 
nant. The integration in ( p^ ) is analytic for the example 
&, giving 

i^HiO ~e + mii/a. (16) 

Paradoxically, vH[£) plays an important role in (jj) at 
1^1 = 0(^/), however small v may be. This is because 
of the contribution from the shoulders of mature shocks. 
If this is so, how can the H term be negligible for larger 
negative ^, as assumed in obtaining (^)? 

Differentiation of (10) gives 



= 3^' + {ul/v)^. (17) 

If Ms = 0{urms)i the sccoud term in (|l7| ) is 0{R^f£,), 
and it dominates the first term for |^| ^ R£,f- Consider 
the measure, per unit length of x, of the portion of a 
shock shoulder where |^| = 0(C/)- If the shock spacing 
is 0(l//c/), this measure is Ms(|C| = O(0)) = 0{R-^) 
while the total measure per unit length Mrd'CI = 0(^/)), 
where |^| = 0{£,f), is 0(1). The total measure is dom- 
inated by field, freshly injected by the forcing, that has 
not interacted with shocks. If R is large, i^i?($ — 
OiS,f ))Q{S, — 0{£,f)) is dominated by the contribution 
(Ms/t^)i^ = 0{R^'f) from the shock shoulders. The corre- 
sponding contribution to vH{S) is 0{R£,'j)Ms/MT- Since 
Ms/Mt = 0{R~^), this implies that the shock shoulders 
make an 0{R^^R^^) contribution to H{£^) for <^ negative 
and 0(C/)- Thus the v term in (H) is the same order as 
the 9(^^Q)/9C term for such ^. It is independent of R as 
R — > oo. 

Next consider the putative range (^). Here the domi- 
nant part of within the shock shoulders is again cx ^, 
while by ( p^ the contribution of the shoulders to Q{£,) 
is cx 1/^. The contribution of the shoulders to H{^)Q{£^) 
is again the dominant one. This contribution thus goes 
like and is independent of to leading order. The v 
term in (^) vanishes to leading order when the ^ differen- 
tiation is performed. This is consistent with the neglect 
of the ly term in deriving (^). 

The regions outside the shoulders have negligible iy£,xx 
at large R. A relative-measure estimate, like that above, 
then implies H{^) cx f for ^/ < < R^^'^^f- H{^) 
grows with |^| throughout the range (^), but neglect of 
d[H{^)Q{^)]/d£, appears to be valid. 
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V. STATISTICAL EQUATIONS FOR VELOCITY 
DIFFERENCES 

Let P{Au, r) be the pdf of u[x + r,t) — u{x, t), where 
r is always taken positive. If z = Au/r, then P{Au, r) — 
Q{z,r)/r where Q is the pdf of z. An equation of motion 
for P is obtained by following Lagrangian trajectories as 
in the derivation ^ of (Q). Two trajectories must be 
followed simultaneously in the present case. The result 
is 



dP 
'dt 
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dAu 



dAv 
dx' 



dx"^ dx''^ I 
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Au 



dP 

dr 



dAu 
dx 



Au)P + F{r) 



d^P 
dAu^ 



(18) 



where F(r) = J^[fix+r,t)-f{x,t)][fix+r,s)-f{x,s)]ds 
and x' ~ X + r. The next-to-last term in (|l8|) expresses 
the change of measure at the points x and x'. It contains 
du{x,t)/dx and du{x' ,t)/dx' . Instead of the gradi- 
ent intensification term in (^, there is now the dP/dr 
term. It expresses the change of the label r carried by 
the velocity difference between the two fluid elements as 
they move closer {Au < 0) or apart (Au > 0). 

An equation for Q equivalent to ( [18[ ) may be written 

as 



dQ 
dt 



= r F{r) 



d^Q 



dz^ 



+ z^^ + zQ- 
dz 



^dQ 

dr 



d[H{z,r)Q] 



dz 



+ M{z,r)Q - V- 

where z = Au/r, 

M{z) ^2z + 2r {dz/dr\z), 

H{z) = {uxx{x' ,t) - Uxx{x,t)\z) jr. 



(19) 



(20) 



In (g^, M has been re-expressed by use of the identity 
du/dx = dAu/dr, where u = ^[u{x, t) + u{x', t)]. M can 
be further tranformed |l^ by use of the homogeneous- 
field identity |l|l 



{dz/dr 



1 



Q{z,r) 



dQ{z', 



'-dz'. 



(21) 



This expresses the advective contributions to ( fi^ ) en- 
tirely in terms of Q, at the expense of creating an integro- 
diffcrcntial equation. An equation equivalent to ( p^ was 
obtained by Polyakov Q by exploiting translation invari- 
ance of the equation of motion for the characteristic func- 
tion (Fourier transform) of Q. 

If the limit r is taken, 2; — > ^ and Q{z, r) Q{Oi 
which is independent of r. Also, H{z,r) H{^) and 



'-F{r) B. Then (|T|) goes into (g). 



If r is greater than typical shock widths and z is large 
enough, the right side of (^ is dominated by contribu- 
tions from shocks such that one of the two points x and 
x' lies in the shock front while the other lies outside. The 
integration limits in (^ij) express the fact that all shocks 
with jumps > Au contribute. 

The pdf of z = Au/r should approximate that of f 
if, in the line-segments of length r that contribute sig- 
nificantly at a given z, ^ fluctuates negligibly over the 
distance r. For ^ > or |^| a plausible sufficient 

condition for this is rkj <C 1. For — ^ ^ the typi- 
cal spatial scale of variation of ^ has been decreased by 
squeezing to ^ O/^/l^l- requirement for negligible 
variation of ^ over r is then r <C which can be 

rewritten as rkf <C O/I^l ''I'CI '"rms- This excludes 
Au values due to the presence of a typical shock within r. 
The condition for the pdfs of z and ^ to fall on each other 
throughout the —3 range is consequently rkf <C R~^/'^. 
If i?r^/2 > rfc/ > the collapse of the pdfs of z and 
^ should extend into the —1 range. In general, if rfcy 1 
the pdf of z should fall below that of ^ at negative ^ large 
enough that rkj <^ is violated. 



VI. COMPARISON WITH SIMULATIONS 



Computer solution of (0) by an SX4 machine was car- 
ried out on a cyclic domain of A*" = 2^^ to N = 2"^^ 
points with unit spacing. The forcing had a wavenum- 
ber spectrum of form (A) fc^ exp[— (fc/fc/)^] or (B) 
k^ exp[— (/c/A:/)^]. The forcing field at each time step was 
an independent realization of Gaussian statistics. The 
initial spectrum for u had the form fc^ exp[— (fc/fc/)^], 
with variance chosen to minimize transients in Urms- The 
simulations were continued until the statistics of interest 
were stationary 10^ time steps). The steady state 
Reynolds numbers R — Urms/vkf ranged from 15 to 
18000. The integrations were performed using second- 
order schemes in x and t. Table I shows the simulation 
parameters for the runs reported here, (i?) is the approx- 
imate time average of R over the period in which Q{S,) 
exhibited a statistically stationary state. Statistics were 
averaged over sets of similar runs that differed only in the 
random numbers used in realizing the forcing and initial 
velocity spectra. For Runs 3-5, additional averaging was 
performed over time in the statistically stationary state. 

^1/2, 



Fig. 1 shows logjoL 

,1/2 



[(f) QiO] plotted against 

for R - 15,1200,18000 (Runs 3-5). Note 
the increase in sharpness of the peak as R increases. 
Fig. 2 shows the central region of £,fQ{£,) plotted against 
^/^/ for three runs with forcing spectrum (B), R ^ 
15, 1200, 18000 (Runs 3-5) and Run 1, with forcing spec- 
trum (A) at i? ~ 15. With the ^/ scaling, the central 



4 



region is substantially insensitive to change of forcing 
spectrum or R. 

Fig. 3 shows superimposed plots, for R ^ 1200 (Run 
4) and R - 18000 (Run 5), of logioK/Q(0] against 
logd^l/^/) for ^ < 0. The straight lines have slopes of 
—3 and —1. The R ~ 18000 data seem consistent with 
the existence of the proposed asymptotic ranges and 
( p^ ) but clearly even higher R values would be needed to 
make an unambiguous case. 

The steady-state pdf equation, obtained by setting 
dQ/dt = in (Q) has interesting stability properties. 
If the dissipation term i/d[H{S^)Q{£^)]/dS, is omitted, the 
general solution is analytically accessible. It contains ex- 
ponential factors that make numerical solution violently 
unstable if carried out in the negative ^ direction and 
highly stable if carried out in the positive ^ direction. A 
consequence is that the central peak of Q is quite insen- 
sitive to the form assigned to the v term. 

Fig. 4 shows the central part oi ^fQ{£,) plotted against 
^/^/ for four cases: (a) the R ~ 1200 simulation (Run 4); 
(b) the left-to-right numerical solution of (||) in steady 
state with dissipation term set to zero; (c) the left-to- 
right numerical solution with dissipation term taken as 
0.45^/(5('C); (d) the left-to-right solution with dissipation 
term taken as 0.8609^/^(0/(1 + CVC/)- The three nu- 
merical solutions were started at large negative ^ where 
QiO ^ and normalized to unit probability. 

The first notable thing about Fig. 4 is that case (b) is 
so close to case (a). It is has a small, unphysical negative 
region at ^ > 0, but nevertheless lies near (a) over the 
peak region. The forms (c) and (d) for the dissipation 
are even closer to (a). They decay according to (|^) as 
^ — > oo. 

Case (c) is the case a = of the closure approximation 
i'd[HQ]/dS, = (a^ + b^f)Q introduced by Polyakov |] 
and studied also by Boldyrev This closure was ob- 
tained as a low-order truncation of the operator-product 
expansion associated with Burgers equation. 

For our present purposes, (c) and (d) are simply two 
functional forms that serve to increase d^Qjd^ for |^| ^ 
i^j. The choice of a functional form for this purpose is 
only weakly constrained, but the numerical coefficient is 
not. It is fixed by the requirements that (3(0 be positive 
for all ^ and decay faster than algebraically as ^ ^ oo 
|^,|2|. The second requirement is satisfied also by case 
(b). If the numerical coefficients in cases (c) and (d) are 
increased, the right tail of Q(0 decays like instead 
of like (|^). If the coefficients are decreased, a negative 
region like that in case (b) is produced. The numbers 
0.45 and 0.8609 stated above are approximations to the 
exact marginal values. 

It was noted in Sec. V that viscous effects from the 
shoulders of mature shocks are present at |^| ~ ^ even as 
R — !■ oo. Fig. 4 suggests that the principal consequence 
of these effects in the Q equation are not upon the form 
of the central peak, whose shape is remarkably stable, 



but upon the decaying region of Q at larger positive ^. 

The putative asymptotic range (|^) is difficult to define 
well by simulation. Q(0 decays so rapidly with increase 
of ^ that very large sample sets are needed. Fig. 5 is a 
plot of 3B9(lnQ)/9(^3) ^ZjQ foj. 2, a i? ~ 15 

simulation with forcing spectrum (B). The exp(— ^'^/3B) 
factor in (^ seems supported. But the approach to the 
horizontal asymptote is protracted, consistent with the 
presence of a positive-power prefactor. 

A definitive test of the prefactor exponent in (^ is 
more difficult. Fig. 6 is a plot of 39(ln Q) / 9(ln ^) against 

jB. If the exp(— ^■^/3i3) factor is present in Q, the pref- 
actor exponent is given by the intercept at ^'^ = of a 
straight line drawn through this plot at large • 

In order to help resolve the prefactor, we have included 
in Figs. 5 and 6 curves corresponding to a mapping ap- 
proximation [|l[ carried out with the same v and forcing 
parameters as the simulation. A detailed description of 
the approximation for the forced case will be given else- 
where. We report now is that it yields (0) at infinite R 
but gives Q(0 oc C/^C exp(— const ^O^rms ~ for 
C ^ finite R. The simulation and mapping curves 
in Fig. 6 lie close to each other. Straight lines drawn 
through the outer parts of both curves (10 < jB < 15) 
intercept the vertical axis near 1. 

It is also difficult in the simulations to resolve the con- 
ditional mean H{^) cleanly for large negative In ad- 
dition to the need for large sample sets, the x grid must 
be sufficiently fine to resolve unusually narrow shocks. 
Fig. 7 shows vH{C)l^) plotted against ^/^ {£, < 0) for 
Run 3, which has high resolution. Also shown are the 
parabolas vH{£,) — and vH{£,) = + Ccf, where £,c 
is the value at which the simulation data for H{S) change 
sign. The latter function is an approximation suggested 
by (|l7|). The second term in the asymptotic relation (|l5| ) 
obviously is negative for ^ < 0. This implies that is 
an upper bound to vH(^) for large negative 

Fig. 8 shows data for ^j:rP{Au, r) from the R - 1200 
simulation (Run 4) plotted against Au/r^f for a num- 
ber of values of r. Also shown is £,fQ{^) plotted against 
£,/£,/■ As r decreases, the curves for P follow that for Q 
over an increasingly long range of Au/r. At the smaller 
r values, this collapse extends into the —1 range of the 
Q curve. The envelope of the knees of the P curves, be- 
yond the region of collapse, follows a line of slope —2 
on the log- log plot. This is because the regions of large 
Au/r are dominated by contributions from the mature 
shocks, and scale with Urms- The latter scaling is clearly 
shown in Fig. 9, where P{Au,r)/£^fr is plotted against 
Au/urms- The outer regions of the curves collapse. 

Figs. 10 and 11 are similar to Figs. 8 and 9 except that 
they show data for R - 18000 (Run 5). 



5 



VII. DISCUSSION 

The theory of Burgers equation with low-wavenumber, 
white-in-time forcing has had a variety of treatments, and 
there has been a variety of predictions. The exp(— ^^/SB) 
factor in the right tail of Q{£,) seems generally accepted. 
It was first predicted by Polyakov 0] on the basis of the 
truncation vd[HQ]/d£_ = {a^ + b^f)Q of the operator 
product expansion, as noted in Sec. VI. Later it was 
recovered by instanton analysis ||^[lOt| . 

Equation (^) is exact. If the i/ term in (Q) can be ne- 
glected in the right tail, it follows that the prefactor is ^, 
as shown in . Neglect of viscous effects in the right tail 
is plausible but not obviously justified. One must con- 
sider the effects of possible proximity of regions of large 
positive ^ to the shoulders of strong shocks. 

All values a 7^ in Polyakov's closure give significant 
viscous effects in the right tail, with the consequence that 
the predicted prefactor is Instanton analysis that 

neglects viscous effects at the start must yield the ^ 
prefactor. However, extraction of the prefactor from in- 
stanton analysis is a delicate matter that requires careful 
treatment of fluctuations about the saddle point. 

There is also disagreement about the exponent of the 
putative power-law range for intermediate negative ^. 
The predictions have rang ed from -2 [|j7| to -7/2 [0. 
Polyakov's closure gives g = 3 — a; values < a < 1 have 
been considered 

As discussed in Sec. IV, the dissipation contribution 
from the shoulders of a distribution of ideal Burgers 
shocks gives iJ(^) oc |^|^ if Q{(,) cx so that the 

1/ term in (^) vanishes to leading order, as required for 
consistency. If Q{£,) oc |^|^'' (q 7^ 3), (||) in steady state 
requires 

^ (22) 
z q 

Thus H{£_) is negative if 2 < g < 3 and positive if g > 3 
OY q < 2. If i? is large, contributions to H for the range 
1^1 ^ R^f in question can come only from the shoul- 
ders of shocks. The profile of ^ in a shock is a negative 
peak so that, whatever the precise profile, the curvature 
of the shoulders is negative and the contribution to H{£^) 
therefore is negative. This appears to rule out q > 3 and 
q < 2. Values q ^ 3 imply shock profiles different from 
that of the ideal Burgers shock (0). 

Values q < 2 are ruled out for another reason. Homo- 
geneity requires J^^QiO^d.^ = 0- The contribution to 
this integral from ^ > is finite and independent of R as 
R ^ 00, if the right tail of Q falls off rapidly in the fash- 
ion generally accepted. But if the negative-^ power-law 



range extends from |^| = 0(^/) to \^\ = 0{R'^^f) where 
c > 0, then the contribution of this range to the inte- 
gral is negative and becomes infinite as i? — > 00. Since 
the contribution from the left tail beyond the powerlaw 
range is also negative, the homogeneity condition cannot 
be satisfied in the limit. 

The value q = 7/2 was proposed in the course of an 
analysis of the inviscid Burgers equation [ pT| . One thing 
established rigorously in this work is that the mean spac- 
ing of strong shocks is 0(l/fc/) if the forcing has compact 
spectral support about kf. The value q — 7/2 was ob- 
tained by examining the behavior of it (a;, t) in the imme- 
diate vicinity of the formation of an incipient shock, and 
then looking back in time at the structure of the regions 
that were destined to form these vicinities. We believe 
that this procedure yields a biased sample of the pre- 
shock velocity field. Equation is intended to describe 
the probability balance associated with the evolution of 
all regions of negative ^. Only a zero-measure set of these 
regions is destined to form part of the immediate vicinity 
of an incipient shock in the limit i? — )■ 00. 

One effect not taken into account in the derivation of 
is passage of a mature shock through a region where 
^ is steepening. Such passage can wipe out the local 
process. Since negative ^ increases in magnitude until it 
either forms part of a shock or is wiped out by a shock col- 
lision, the collision process should make Q{£,) fall off more 
rapidly than (H) with increase of |^|. The rate of increase 
of within a fiuid element, before viscous effects are felt 
and neglecting forcing effects, is £,{t) oc ^o/(l — \^o\t), 
where t is measured from some effective initial time when 
^ = ^0- The times at which |^| is large compared to |^o| 
are all crowded into an interval <C 1/ICol at i — O(l/|^o|). 
The mean time between shock collisions is 0{1/ kfUrms)- 
Because of the crowding in time, the the collision prob- 
ability changes very slowly with |^| for |^| ^ and we 
do not expect the collision effects to change the power 
law in (||). 

We are indebted to S. A. Boldyrev, S. Chen, A. 
Chekhlov, M. Chertkov, G. Falkovich, V. Lebedev, A. 
M. Polyakov, B. Shraiman, Y. Sinai, and V. Yakhot for 
valuable discussions and correspondence. T.G.'s work 
was supported in part by Grant-in- Aid for Scientific Re- 
search C, 09640260 from the Ministry of Education, Sci- 
ence, Sports and Culture of Japan, and by the Computa- 
tional Science and Engineering, Research for the Future 
Program, Japan Society for the Promotion of Science. 
R.H.K.'s work was supported by the United States De- 
partment of Energy, Office of Basic Energy Sciences, un- 
der Grant DE-FG03-90ER14118. 
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TABLE I. SIMULATION PARAMETERS 



Run 


(R) 


Urms — 0) 


Spectrum 


kf 


V 


iV 


dt 


B 


Set size 


1 


15 


1 


A 


0.02 


10/3 




0.1 


4 X 10"" 


102 


2 


15 


1 


B 


0.02 


10/3 


2" 


0.1 


7.2 X 10"® 


102 


3 


15 


1 


B 


0.02 


10/3 


220 


0.1 


7.2 X 10"® 


102 


4 


1200 


1 


B 


5 X 10"* 


2 


2l8 


0.1 


2 X 10~i° 


102 


5 


18(KK) 


J. 


B 


5 X 10 


2.2 


2-" 


0.2 


X 10 " 


()1 
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FIGURE CAPTIONS 

FIG. 1. Plot of log]^Q[(^^)^^^ Q(^)] against ^/ (^^)^^^ for three sets of runs with forcing spectrum (B) and R ^ 15 
(Run3, hght dashes), 1200 (Run 4, heavy dashes) and 18000 (Run 5, solid line). 

FIG. 2. The central region of C/Q(C) plotted against for the three runs with forcing spectrum (B), R ^ 15 (Run 
3, light solid line), 1200 (Run 4, heavy solid line), 18000 (Run 5, heavy dashes), and a run with forcing spectrum (A) 
at i? 15 (Run 1, light dashes). 

FIG. 3. Plot of logio[C/Q(0] against log(|C|/^/) for ^ < 0. i? - 15 (Run 3, light dashes) R - 1200 (Run 4, heavy 
dashes) and R ^ 18000 (Run 5, solid line) are shown. The straight lines have slopes of —3 and —1. 

FIG. 4. The central part of £,fQ{£,) plotted against ^/C/ for four cases: (a) the R ^ 1200 simulation (Run 4, heavy solid 
line); (b) the left-to-right numerical solution of (|^) in steady state with dissipation term set to zero (light solid line); 
(c) the left-to-right numerical solution with dissipation term taken as 0.45^/(5(^) (heavy dashes); (d) the left-to-right 
solution with dissipation term taken as 0. 8609^/(5(0/(1 + (.'^ / ^i'}) (light dashes). 

FIG. 5. Plot of 3Bd(lnQ) /d{£,^) vs £,^/B for a ~ 15 simulation (Run 2, solid line). Dashed line is a mapping 
approximation. 

FIG. 6. Plot of 39(lnQ)/9(ln^'^) against / B (Run 2, solid line). Dashed line is a mapping approximation. 

FIG. 7. Plot of vH{£)/^'j against negative values of ^/^f for a i? ^ 15 simulation (Run 3, solid line). Also shown are 
the parabolas {(./^f)"^ (light dashes) and (^^ +^cO/?/ (heavy dashes). 

FIG. 8. Plot of ^frP{Au, r) from the R ^ 1200 simulation (Run 4) against Au/r^f (dashed lines) for values of r that 
increase by factors of \/2 from rkf fa .008 to rfc/ « 0.72. The curve that is highest on the left side is rfcy « 0.72. Also 
shown is ^fQ{£,) plotted against ■^/C/ (solid line). 

FIG. 9. P{Au, r)/£^fr plotted against ud/urms for the R ^ 1200 run. The curves denote the r values of Fig. 8. 

FIG. 10. Plot of ^frP{Au,r) from the R ~ 18000 simulation (Run 5) against Au/r^f (dashed lines) for values of r 
that increase from rkf k, 0.0016 to rkf k, 0.82. The curve that is highest on the left side is rkf « 0.82. Also shown 
is ^fQiO plotted against (soHd line). 

FIG. 11. P{Au,r) I ^jr plotted against Au/urms for the R ^ 18000 run. The curves denote the r values of Fig. 10. 
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